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Abstract 

A nonlocal Landau-type free energy functional of the spin density is developed to model the large 
variety of magnetic states which occur in the magnetic field-temperature phase diagram of magne- 
toelectric CuFe02- Competition among long-range quadratic exchange, biquadratic anti-symmetric 
exchange, and trigonal anisotropy terms, consistent with the high-temperature rhombohedral R3m 
crystal symmetry, are shown to all play important roles in stabilizing the unusual combination of 
commensurate and incommensurate spin structures in this highly frustrated triangular antiferro- 
magnet. It is argued that strong magnetoelastic coupling is largely responsible for the nonlocal 
nature of the free energy. A key feature of the analysis is that an electric polarization is induced by 
a canting of the non-collinear incommensurate spin structure. Application of the model to ordered 
spin states in the triangular antiferromagnets MnBr2 and NaFeO"2 is also discussed. 
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I. INTRODUCTION 



One of the more intriguing features of the magnetic field-temperature (H — T) phase dia- 
gram of the rhombohedrally stacked triangular antiferromagnetic CuFe02 is the occurrence 
of a field-induced non-collinear incommensurate (IC) spin structure amongst a multitude of 
commensurate collinear magnetic ordered states.— ' 2 » 3 ' 4 ' 5 ' 6 ' Usual Hund's rules suggest that the 
magnetic Fe 3+ ions in this semiconductor do not exhibit low-order spin-orbit coupling since 
L = (with S = 5/2). The origin of the observed c axis magnetic anisotropy is not precisely 
understood but it is generally agreed that it should be very weak.->& The frustration inherent 
in triangular antiferromagnets is known to give rise to period-3 (P3) elliptically or helically 
polarized (S) spin structures at H = 0. A field applied along the anisotropy (c) direction 
typically induces a spin flop transition S _L c in weak axial antiferromagnets.-^ 1 ^ 1 ^ In the 
case of CuFe02, spin non-collinearity occurs only at moderate values of H, and high-field 
spin states are linearly polarized with S||c. Key to the understanding of these unusual fea- 
tures is the very high degree of frustration resulting from not only the triangular geometry 
but also the long-range exchange interactions, magnetoelastic coupling, and anti-symmetric 
exchange leading to the magnetoelectric effect. In our previous work (Ref. 
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hereafter 



referred to as I), the observation of the field-induced sequence of period-4 (P4), IC, period-5 
(P5) and P3 basal-plane modulated spin structures at T = was shown to be a consequence 
of competition among these disparate interactions. As the temperature is lowered in zero ap- 
plied field, CuFe02 exhibits successive magnetic transitions at T/vi ~ 1AK to a collinear IC 
state followed by a discontinuous transition at T/y 2 — UK to a collinear P4 spin structure. 
A similar sequence of transitions has also been reported in the triangular antiferromagnets 
MnBr 2 and NaFe0 2 .^^ 

In the present work, a representation of these effects is incorporated into a non-local 
Landau-type free energy functional of the spin density, F[s(r)], that is constructed from 
symmetry arguments, in order to develop a model of the complex H—T phase diagram.- The 
formalism is essentially phenomenological but contains the same type of T = interactions 
that are considered in usual spin Hamiltonians. In the case of CuFe0 2 , these include in- 
plane exchange interactions up to third neighbor, Ji, J 2 and J^&H^ML as well as inter-plane 
exchange J'.2J^21 These interactions lead to a minimization of the wave vector dependent 
exchange interaction Jq near a multicritical point where multiple periodicities are close 
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in energy. In addition to various biquadratic exchange interactions (symmetric and anti- 
symmetric) and usual axial anisotropy, the rhombohedral R3m crystal symmetry also allows 
for the existence of an unusual trigonal interaction which was previously used to explain 
the magnetic structure anomalies in pure Ho.™ The importance of this term in stabilizing 
canted structures is explored here. This formalism also leads to umklapp-type terms in the 
free energy which are non-zero only if the ordering wavevector is equal to a multiple of the 
crystal reciprocal lattice vector G, i.e., nQ = G with n having values 3, 4 and 5 in the 
present case. Within the Landau formalism, these types of terms are then responsible for 
the stability of the commensurate phases depending on the values of T and H . A Landau- 
type expansion based on molecular field theory applied to MnBr 2 is presented in RefJl4 



which contains features in common with the present model. A more formal group theoretic 
approach to understanding spin structures in a variety of multiferroic compounds, with a 
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discussion of CuFe02, is given in Ref. 

The nonlocal formalism naturally leads to fourth-order, and higher, contributions (in s(r)) 
to the free energy which are dependent on Q and illustrates explicitly that terms which are 
independently invariant (with respect to all symmetries) have independent coefficients.— ^ 
This contrasts with the usual local form of the Landau energy which results in all isotropic 
terms having the same coefficient at a given power of s. From a phenomenological point of 
view, the strong symmetry arguments are sufficient to assign different values to each of the 
many coefficients. However, such nonlocal effects can also be understood to have microscopic 
origins by considering interactions between magnetic and other degrees of freedom, such 
as magnetoelastic coupling, which gives rise to nonlocal biquadratic exchange terms. Spin- 
lattice coupling is known to be strong in CuFe02,— <2£i2&i2£i which provides further justification 
for the present approach. 

II. NONLOCAL LANDAU FREE ENERGY 

The development here of a nonlocal Landau-type free energy functional of the spin density 
F[s(r)] for CuFeC>2 follows the description given in Ref. llO for triangular antiferromagnets 
and uses the symmetry arguments given in I. Although the expansion is carried out to sixth 
order, a number of simplifications can be made for the purpose of understanding the phase 
diagram. It is convenient to write our model for F as the sum of isotropic and anisotropic 
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contributions in the form 



F = F 2 + F 4 + F 6 + F z + F CP + F K - m • H (1) 
where the first three terms are the isotropic contributions relevant for all magnetic systems 

F2 = w*J drdr ' A ( r - r >( r ) ■ s ( r ') ( 2 ) 

F 4 = J dridr 2 dr3dr4-B(ri,r2,r3,r 4 )[s(ri) • s(r 2 )][s(r 3 ) • s(r 4 )] (3) 
F 6 = J dn • • • dr 6 C(n, • • r 6 ) [s(r x ) • s(r 2 )] [s(r 3 ) • s(r 4 )] [s(r 8 ) • s(r„)] . (4) 

Temperature dependence enters in the usual Landau treatment of the isotropic second-order 
term, here generalized to be of the form 

A(r) = ak B T5(r)+fj(r) (5) 



where a depends on the total angular momentum number j and J(r) is the exchange inter- 
action. Such a form can be derived from a molecular field treatment of the corresponding 
Heisenberg hamiltonian.— i2L2£ This mean field treatment also yields local forms for the higher 
order isotropic terms, involving constants B and C which also depend on j and are propor- 
tional to k B T. The nonlocal form of spin interactions, such as biquadratic exchange, can 
arise from a variety of n-body interactions^ or indirectly from the coupling of s(r) to other 
relevant degrees of freedom such as lattice or electronic (see the Fqp term below). Nonlocal 
fourth order contributions to the energy of the form given above that arise from magnetoe- 
lastic coupling have been derived.— 1 ^ In the present model, we set k B = 1, j = 1 and simply 
treat a and the nonlocal fourth- and sixth-order coefficients as phenomenological parameters 
of undefined microscopic origin but with the support for their existence in CuFeC>2 from the 
strong spin-lattice coupling and a nontrival spin-polarized electronic structure.— 

There are a large number of independent anisotropic contributions invariant with respect 
to the R3m crystal rhombohedral symmetry group generators {5*61000} and {<t„|000}.— Only 
the essential terms are considered here. Axial anisotropy (z||c) is included only at second 
order in s: 

Fz = W~J drdr ' J ^ r - r 'K(rK(r'). (6) 

Note that although this general form includes both single-ion anisotropy D, where r = r', 
as well as two-site anisotropic exchange, the model used in I which includes only the latter 
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is adopted here. Since the anisotropy is small, such a distinction has little impact on the 
results presented below. Higher order terms such as (s ■ s)sl and are omitted here for 
simplicity. 

The stability of the field-induced non-collinear phase at zero temperature was shown in 
I to be a result of including a biquadratic anti-symmetric exchange term of the form 

Fcp = -^X~ p E / drdr ' [C^HsM * s(r')) ■ zf (7) 

where r = r — r' and a = x,y. Such a term is invariant with respect to all relevant 

symmetries and therefore must exist in the energy. Contributions to it can also 

result of coupling between the electric polarization, lattice and spin vectors in the forrn^ 

Fc = W^S drdr ' c ( r )l p ( T ) x f ] ■ x < r % ( 8 ) 

in combination with the lowest-order contribution to F in P given by^^ Fp = 
-^Ap J drdv' P(t) 2 . Such cross product interactions are known to stabilize helically po- 
larized spin structures.— Further analysis of the dependence of P on s(r) is given in the 
following section. 

Systems with trigonal symmetry also support an unusual anisotropy term which couples 
in-plane and out-of-plane components of the spin vector.— For simplicity we use the local, 
single-site expression of this effect with a coefficient K: 

Fk = wI ^( r K( r )[ 3 *) - *J(r)]. (9) 

This interaction favors linear spin configurations which are canted, having both z and basal- 
plane components. The combination of Fcp and Fk contributions to the free energy thus 
favor a canted elliptically polarized state which induces an electric polarization P in the 
basal plane. 

III. WAVE VECTOR REPRESENTATION 

Positions within the rhombohedral lattice of the Fe 3+ magnetic ions are described here 
through an equivalent simple hexagonal unit cell with three triangular layers A, B, and C, 
using basis vectors wa = 0, w# = |ax + |6y + |cz, wc = \ak — |6y — |cz, with a and 
c being the lattice constants, b = (v3/2)a and x _L y _L z. As discussed previously,— 
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long-range magnetic order is characterized by the function p(r) through the spin density 
written as 

s(r) = ^EE/^(r)<Kr-R) (10) 

where R = Ri + Wj, Ri = (n + |m)ax + mby (with n and m being integers) specifies 
hexagonal Bravais lattice sites, j = A, B, C, and N = 3Ni is the number of magnetic ions. 
The variety of spin structures in the magnetic phase diagram of CuFeC>3 can be adequately 
distinguished with p expressed as a Fourier expansion with a small number of terms in the 
form 

Pj (r) = m + S je iQr + S*e- lQ - r (11) 

where m is the uniform magnetization induced by an applied field along the z axis, S is the 
spin polarization vector and Q is the primary wave vector. A full description of some of the 
ordered states requires additional Fourier components, as discussed below. 

With only near-neighbor coupling J' between planes, a reasonable ansatz is to assume 
that the complex spin polarization vectors Sj on adjacent layers differ only by a simple phase 
factor (p with each layer having an overall phase factor 7 in the forrn^ 

S A = Se 17 , S B = S C = Se i(7 ~^. (12) 

The primary description of the spin polarization of each state is further characterized by 
writing^ 

S = S!+*S 2 (13) 

where Si and S2 are real vectors. Previous analysis of Landau-type models to describe the 
magnetic states of frustrated triangular antiferromagnets&^ii 1 ^ have demonstrated that 
sufficient flexibility can be achieved by assuming Si and S2 are characterized by writing 

Si = Scos{3[sm9pi + cos#z], S 2 = 5 , sin/5p 2 (14) 

where pi _L _L z, i.e., p\ and p2 are in the hexagonal basal plane. The Q — th component 
of the spin density on layer A is then given by 

s A (r) = 2Si cos(Q ■ r + 7) - 2S 2 sin(Q • r + 7). (15) 
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Linearly polarized states are described with S2 = 0, proper helically polarized spin config- 
urations are realized by Si _L S2 J- Q with Si — S2 (f3 — 7r/4), and elliptically polarized 
states have Si S^ ^ 0. It is also convenient to define dimensionless wave vectors in units 
of the lattice constants 

q x = aQ x , q y = bQ y , q z = cQ z . (16) 



The zero-field till P4 phase described in Fig. 11 of Ref. 
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for a single layer is characterized 



by the above relations with q x = n, q y = 0, or equivalently q x = tt/2, q y = 37r/4 due to 
the triangular symmetry, along with j3 = 0, 9 = and 7 = 7r/4 (also see Refill). In terms 
of triangular lattice basis vectors a = ax and b = |ax + by, the relations are q a = q x and 
Qb = \lx + q y - Thus, for example, the above degenerate modulations are also described by 
Qa = 7r, q b = tt/2, and q a = tt/2, q b = it. 



A. Second Order Isotropic Terms 

The free energy density as a function of m, S and Q is developed by using the spin 
density given by (10) and (11) in expressions (l)-(9). As shown previously—^ this analysis 
naturally leads to umklapp terms arising from the condition for periodic lattices: 

^£e inQRl =A nQ>G (17) 

where G is a hexagonal reciprocal lattice vector. Second order isotropic contributions reduce 
to 

F 2 = \A Q m 2 + AqS 2 (18) 
where S 2 = S • S*, Aq = aT + Jq and 

^Q = ^E<W QR (19) 

iV R 

Within the present model, this leads to the following wave vector dependence of the exchange 
integral Jq = 2/(q, 0), with 

/(q,0) = JMq) + J 2 / 2 (q) + J 3 / 3 (q) + ±J'f '(q)(l + 2cos0) (20) 
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whereas 

fx = cos q x + 2 cos \q x cos q y 

f 2 = cos 2g y + 2 cos |g x cosg^ 

/ 3 = cos 2^2; + 2 cos g x cos 2g y 

/' = cos(|g x - \q z ) + 2 cos |g x . cos (\q y + |g 2 ) (21) 

and J\ = J (ax), J 2 = J(26y), J3 = J(2ax), J' = J(w). Note that Jj > represents 
antiferromagnetic coupling. Using these results, the coefficient of m 2 can be expressed as 
A = aT + 2/(0,0). 

Within this mean field theory, the wave vector characterizing the first ordered phase to 
appear as the temperature is lowered from the paramagnetic state (S = 0) is determined 
by the extrema of Jq (the effects of anisotropic exchange are discussed below). Results of a 
numerical algorithm to sketch the J 2 — J3 phase diagram are shown in Fig. 1 (also see Ref. 



36l ). Here, we set J\ = 1 for convenience and consider two values of inter- layer coupling, 
J' = and J' = 0.4, as in I. In the case of J' = (broken lines), the usual P3 modulation 
found in the frustrated triangular antiferromagnet (associated with, for example, 120° spin 
structures) occurs in the upper left part of the diagram. At more positive values of J 2 , 
further frustration stabilizes an IC modulation. Simple antiferromagentic structures^ (P2) 
are found in the lower right region. Note that due to the triangular symmetry, there are a 
number of equivalent wave vector descriptions for the same basic structure. For example, 
120° spin configurations with (47r/3,0), (27r/3,7r) differ only in their chirality.— Period- 
2 structures with (0, 7r), (27r, 0) and (tt,tt/2) differ only by a rotation of axes. With the 
addition of interlayer exchange J' = 0.4 (solid lines), there is a nontrivial interplay between 
q z and q x , q y as can be seen in Eq. (21). The boundary between P3 and IC phases disappears 
and there is an additional AF or IC modulation between adjacent planes. For example, in 
the lower part of the phase diagram, degenerate modulations (0, 7r, 2tt) and (2tt, 0, 0) both 
yield simple AF inter-plane structures since q • = it and q • wp = it. Also indicated on 
Fig. 1 is the point corresponding to J 2 = 0.3, J3 = 0.3 (close to the values in I) which are 
used in the calculation of the H — T phase diagram described below. 
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FIG. 1: Sketch of the J 2 — J3 phase diagram based on minimization of the exchange integral J q 
with J\ = 1. Broken curves correspond to the case J' = and solid curves J' = 0.4. Solid circle 
indicates values used in the present model: J2 = J3 = 0.3, J' = 0.4. 

B. Fourth and Sixth Order Isotropic Terms 

Higher order terms can also be evaluated with the assumed spin density (10)-(12). For 
convenience, regular (non-umklapp) terms and umklapp terms are written separately. In 
the case of the isotropic fourth-order contributions, the result can be expressed as follows: 

F4 = F^r + F 4j3 + F 4i4 (22) 

where 

F AA = BiS 4 + |S 2 |S • S| 2 + \B 3 m 4 + 2B 4 \m ■ S| 2 + B 5 m 2 S 2 (23) 

F 4i3 = S 4>3 [(m • S)(S • S)e 3 ^ + c.c.]A 3Q , G (24) 

F 4 , 4 = |S 4 , 4 [(S • S)V^ + c.c.]A 4Q)G . (25) 
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As shown in Ref. 
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umklapp terms involving 2Q are accounted for within the regular terms 
by a suitable renormalization of the spin density amplitude. Expressions for some of the 
coefficients in terms of the Fourier transform of the nonlocal function 2?(r 1; r 2 , r 3 , r 4 ) are 



given in Ref. 
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and the others may be easily deduced. For example, Bi = £?q,-q,q,-q, B a = 
-£>o,q,o,-q and -B43 = -£>o,q,q,q- The important point here is that the nonlocal formulation 
naturally leads to the result that each independently invariant term has its own independent 
coefficient. In a local formulation, all isotropic terms have equal coefficients. In the present 
model, each of these fourth-order (and sixth-order) independent coefficients is assumed to be 
constant. As discussed previously,-^ non-collinear spin structures are stabilized with B 2 > 
and collinear states with B 2 < 0. Usual spin-flop transitions in antiferromagnets occur as a 
consequence of having B± > so that a spin configuration with S _L H is stabilized. 

The above umklapp terms are non-zero only if 3Q = G (and m 7^ 0) or 4Q = G, 
allowing for the possibility of energy reduction if the system assumes these periodicities.— 
In the case of S||m for the first term, or in the case of any collinear state for the second 
term, these expressions reduce to 2£> 4j3 mS' 3 cos(37) and ^B^S 4 cos(47), respectively. For 
positive coefficients, these two terms are then each minimized by phase factors 7 = 7r/3 and 
7 = 7r/4, respectively. 

In an effort to reduce the number of model parameters, a somewhat simplified approach 
is adopted to treat the isotropic sixth order contributions. In addition to the regular terms, 
there are umklapp terms involving periodicities G/3, G/4, G/5 and G/6 . For simplicity, 
it is assumed that each of the independent terms forming the regular part have the same 
coefficient C. This is equivalent to a local formulation of these contributions. Umklapp 
terms, however, are each assigned an independent coefficient. The result can be expressed 
as 



^6 — F&,R + -^6,3 + ^6,4 + ^6,5 + -^6,6 (26) 

where 

F^ R = ±C{S$, + 6S*\S ■ S| 2 + 24S 2 |m • S| 2 + 12[(S ■ S)(m • S*) 2 + c.c.]} (27) 

^6,3 = |C 6 , 3 {[12^(m-S)(S-S) + 6(m-S*)(S-S) 2 + 8(m-S) 3 ]e 3 ^ + c.c.}A3Q, G (28) 

^6,4 = |C M {[3S' 2 (S-S) 2 + 20(m-S) 2 (S-S)]e 4 ^ + c.c.}A 4QiG . (29) 
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(30) 



(31) 



Note that odd-order umklapp terms occur only in the presence of a magnetic field where 
m 7^ and have a maximum effect on lowering the free energy with linearly polarized 
spin configurations and S||H. In-plane and c-axis periodicities, n and m, respectively, can 
be separated by writing Q = ^Gj_ + ^G| where Gj_ and G|| are reciprocal lattice vectors 
perpendicular and parallel to the hexagonal c axis. At low values of magnetic field, antiferro- 
magnetic inter-plane exchange J' > dominates and period-2 inter-layer c axis modulations 
are stabilized. For the odd-period in-plane commensurate structures, Q| = 0, G|| and the 
umklapp conditions 3Q = G and 5Q = G are satisfied for P3 and P5 spin configurations, 
respectively. These characterizations are consistent with neutron diffraction data.->2 

A complete description of the spin structures usually requires that additional wave vector 
components (Q') be added to the Fourier expansion of the spin density (11) which lead to 
a 'squaring-up' of magnetic structures.— 1 ^ Fourth-order isotropic umklapp terms of the 
form (S' • S)(S • S) exist for cases where Q' + 3Q = G leading to the incipient relation 
S' ~ S(S • S). With a magnetic field present, terms of the form (S' • m)(S • S) occur at fourth 
order if Q' + 2Q = G giving S' ~ m(S • S). A larger number of possible secondary wave 
vectors arise from sixth order umklapp terms. A fully consistent analysis would lead to many 
additional contributions to the free energy but are not required for the semi-quantitative 
description of the phase diagram given here. 

C. Anisotropic Terms 

Evaluating the anisotropic contributions to the free energy, F z , Fcp and Fr, in terms 
of the spin density parameters m, S and Q follows the method as described above. Since 
anisotropy is known to be small in CuFeC>2, anisotropic umklapp terms are omitted for sim- 
plicity. The resulting anisotropic exchange terms are identical in form to isotropic exchange, 
but involve only z components of the spin vectors: 



Here, J z q is given by the relations (19)-(21) but with the isotropic exchange parameters Jj 
and J' replaced by their anisotropic counterparts J zi and J' z , as in I. The coefficient J z0 is 



F = - 

L z o 



\JzQ m2 z + Jzq\S z 



(32) 
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then given by this expression evaluated at Q = and (b = in (21). 

An expression for the nonlocal biquadratic antisymmetric exchange interaction can be 
determined by first evaluating Fc using (10) and (11) with the assumption that the electric 
polarization vector P is uniform. This gives^ 

F c = i(C x P x + C y P y )i ■ (S x S*) (33) 

where, as in I, magnetoelectric-type interactions C\ = C(ox) and C = C(w) are included, 
giving 

C x = -|6{Ci cos|g x sinq y - |C"[sin(§<^ - \q y ) - sin^^ + \q y ) cos|gJ(l + 2cos0)} 
C y = |a{Ci(sing x . + sin ^q x cosq y ) + C sin \q x cos(|gj / + |g z )(l + 2 cos0)}. (34) 

Minimization of the free energy F C p = F c + F P , where F P = ^Ap(P% + Py), thus yields the 
relations between P and the spin polarization vectors: 

P x = -{i/A P )C x {$ x S*), 

P y = -(i/A P )C y (S x S*),. (35) 

The combined antisymmetric biquadratic exchange contribution to the free energy then 
takes the form 

FcP = ^- p (C 2 x + C*)(SxSr z . (36) 

These relations make clear that a uniform electric polarization cannot be induced by a 
collinear spin structure where (3 = or 7r/2 since 

(S x S*) z = 2i(Si x S 2 ), = iS 2 sin 2(3 sin 6. (37) 

Note also that P = in cases where S lies strictly in the z — p 2 plane (6 = 0), as has been 
recently emphasized.- 1 ^ 1 ^ 1 ^ 

Finally, the local formulation of the trigonal anisotropy term given above can be expressed 
in terms of the spin polarization vector components as 

F K = K{[3(S* x ) 2 S y S z - S z S y (S;) 2 + 2S y S;(3\S x \ 2 - \S y \ 2 )] + c.c.}. (38) 

Using the parameterizations of S given above, this interaction term can be written as 

F K = 3KS 4 cos 2 (3 sin 2fl(sin 2 (3 - cos 2 (3 sin 2 6) (39) 

which generally favors canted spin structures < 9 < n/2. 
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IV. MAGNETIC PHASE DIAGRAM OF CUFEQ 2 



The nonlocal free energy functional formulated above can be expressed as a function of 
the parameters which characterize the spin density F = .F(Q, S, m, 0, 7, 9, 0). Equilibrium 
spin configurations as a function of temperature and magnetic field are then determined 
by minimization (numerically) of F for a given set of coefficients. Since the phase diagram 
involves spin structures characterized by a number of different wave vectors Q which can 
each have associated umklapp terms, it is necessary to minimize (with respect to the other 
parameters) separately and compare the distinct free energies associated with IC, 3Q = 
G, 4Q = G, and 5Q = G (labeled F IC , F P3 , Fp 4 , F P5 ) phases in order to determine the 
equilibrium state. 

A feature of the present model is the large number of coefficients which is a consequence of 
not only the nonlocal formulation but also the many competing interactions. This multitude 
of effects are, however, essential for a complete understanding of the complex magnetic 
phase diagram of this highly frustrated compound. In order to reduce the number of free 
parameters, zero-T coefficient values from I are used here and other coefficients are a priori 
assigned reasonable values. Only a relatively small number of coefficients are adjusted in an 
effort to reproduce the essential features of the phase diagram. As in I, the overall energy 
scale is set by assigning J\ = 1 with other exchange parameters given by J2 = J3 = 0.3 
and J' = 0.4. Magnetoelectric coupling coefficients used in I are adopted here, G\ = 0.3, 
C = 0.1 and A P — 1, as well as the assignment of an anisotropy strength of 3% so that 
Ji Z = 0.03Jj and J' z = 0.03J'. This level of anisotropy is also adopted here for the trigonal 
coefficient, K = 0.03. Some of the Landau coefficients were arbitrarily set as follows: a — 1, 
Bq — 1, B\ — 1, E>2 = 0.1, -B3 = 0.1 and C = 0.1. Assuming a positive value for the 
term B^\m- S| 2 was always found to yield a low-field spin- flop transition (since anisotropy is 
weak), which is not observed in the magnetic phase diagram of CuFe02 (at least moderate 
field values). Assigning the negataive value B^ = —0.2 serves to enhance the stability of the 
reported configurations with S||H in the P5 and P3 states (also see Ref ill). (This point is 
discussed further below.) Coefficients of the umklapp terms, -84,3, -84,4, C^q = 6*6,3, ^6,4, 
and 6*6,5 are then adjusted to reproduce essential features of the H — T phase diagram. Note 
that the 6Q = G umklapp term in (31) contributes to the P3 state and its coefficient is 
assigned the same value as the sixth-order 3Q = G umklapp term for simplicity. 
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For the range of coefficients considered here, the inter-layer phase factor <fi is always found 
to be zero. The overall spin polarization phase factor 7 appears only in umklapp terms and 
the free energy is minimized with 7 = n/n for commensurate states described by nQ = G. 

Consider first the sequence of transitions which occur at H = 0. With the addition of 
axial anisotropy, the spin structure to first appear as the temperature is lowered from the 
paramagnetic state is linearly polarized S||c. The free energy up to fourth order is given by 

F = (A Q + ^q)^ 2 + {B 1 + \B 2 )S A + \B iA S A cos(4 7 ) A 4Q , G (40) 

with the transition temperature given by T/vi = —(Jq + Jzq)/ '- The wave vector is thus 
determined by the value which maximizes this function, as shown in Fig. 1. For the set 
of exchange parameters used here, this gives the collinear IC phase with degenerate wave 
vectors (q x ,q y ,q z ) = (2<7, 27r, 7r), (2tt — q,3/2q — n,n) and {2ix — 2q, n,7i) (or equivalently 
(q a , qb, q c ) = (2q, 2ir + q, n), {2n — q, q, it), and {2n — 2q ) 2n — q, ir) respectively) where q ~ 0.22, 
giving T/vi — 2.0. The result | < q < | is consistent with neutron diffraction data 2 and 
the speculation of multi-q domain structures.— Numerical minimization of the free energies 
corresponding to IC and P4 (tt, 2tt, tt) phases with all terms included is then performed in 
order to compare Fjc and Fp^ as a function of temperature. 1 - 4 Assigning the parameters 
values -B 4> 4 = 0.8 and C 6)4 = 0.2 is found to yield the result that F Pi < F IC for T < T N2 — 
1.5, which approximately agrees with the experimental data for T^ 2 /Tni- The transition at 
T/vi is continuous while the IC-P4 transition at T/v 2 is discontinuous. 

Using this method of comparing free energies associated with IC and commensurate 
spin states, the phase diagram with H||c is determined numerically. The three remaining 
coefficients are fit to best reproduce the transition boundaries in an effort to achieve semi- 

n 

quantitative agreement with the corresponding experimental results presented in Refsj3| and 
J. This procedure yields B i 3 = 1.0, C 6 6 = 0.1, and C 6i5 = 0.6. The linear IC, P5 and P3 
phases remains stable with S||c at all field values considered due to setting B4 < 0. The 
elliptical IC phase is stabilized at moderate field values with (3 ~ 0.197T at lower T. There is 
little change in wave vector from the zero field values since C 2 and Cy in (36) are small. Due 
to the small trigonal interaction term F Kl this structure is also found to be canted away from 
the c-axis by about 9 ~ 10°. Similarly, at low temperatures, the linear P4 phase exhibits a 
discontinuous transition to canting away from the c axis by about 45° (indicated by the thin 
solid line in Fig. 2). The transition between elliptical and linear IC phases is continuous 
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with (3 acting as an order parameter. The canting of magnetic structures does not occur 
within the present model in the case K — 0. Critical fields at low temperature H P 4_jc — 1, 
Hic-P5 — 2, Hp 5 ^p 3 ~ 3 are in fair agreement with the experimentally observed ratios 
H IC -p 5 /H P4 _ ic ~ 13T/7T and H P ^ P3 /H P4 „ IC ~ 20T/7T.M 




T 



FIG. 2: Sketch of the model phase diagram based for CuFeO"2 on minimization of the free energy (in 
units of Ji). Solid and broken lines represent discontinuous and continuous transitions, respectively. 

V. DISCUSSION AND CONCLUSIONS 

In this work, a model nonlocal Landau-type free energy has been developed and ana- 
lyzed that captures essential features of the complex magnetic phase diagram of the highly 
frustrated magnetoelectric triangular antiferromagnet CuFe02 and serves as an extension 
of a previous zero temperature model.— In contrast with local formulations, the present 
approach naturally leads to the result that each symmetry-invariant isotropic term has its 
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own independent coefficient. In the case of CuFeC>2, it is argued that strong magnetoelastic 
coupling is one source for such nonlocal biquadratic (and higher order) spin interactions. As 
emphasized in Ref. [Lj, magnetoelectric coupling also provides a microscopic mechanism for 
the existence of biquadratic antisymmetric exchange. In addition to weak axial anisotropy, 
the rhombohedral crystal symmetry also supports a somewhat unusual trigonal term which 
couples basal-plane and axial spin components. With analysis performed in terms of a 
Fourier expansion of the spin density, this approach also facilitates comparison with neutron 
diffraction data. 

There are a number of limitations associated with this method. The Landau model is 
an expansion of the free energy around s(r) = so that results far below T/v and at high 
field values may require a large number of terms. In addition, only a small number of mod- 
ulated spin structures can be fully described by just a few Fourier components. The model 
in its present form fails, for example, to account for the very high field spin-flop P3 phase 
suggested by magnetization measurements.- The merging of P4, IC and P5 transition bound- 
aries suggested by the experimentally determined phase diagram is also not found in our 
truncated model. In addition, it is a mean field theory. Within the present formulation, the 
paramagnetic-IC (para-IC) transition temperature is given by Tjvi = — (Jq + Jzq)] 2 /(o&b) 
where a = 3/ [2j (j + 1)]. Using j = 5/2 and the estimate J\ ~ 5.3K(ks) from spin- wave 
data, 20 (along with the other exchange interactions values given above) gives T/vi = 38.6K. 
The nearly factor of 3 discrepancy with the experimental value can be attributed to strong 
critical fluctuation effects associated with the high degree of frustration in this compound. 

The nature of the symmetry breaking at the two continuous transitions, para-IC and 
IC(linear)-IC (elliptical), can be analyzed in a manner similar to the weak axial triangular 
antiferromagnet CsNiCl3. 40 Both should belong to the standard XY universality class involv- 
ing a two component continuous order parameter. In the case of the para-IC transition, the 
two components can be identified as the magnitude of the spin order, S, and the associated 
phase angle 7. In the case of unfrustrated systems, this transition would be identified with 
Ising universality. In the finite-field transition to the elliptical phase, the system develops 
(XY) basal-plane components p\ and f>2- 

A consequence of the present model and analysis is that basal-pane components of the 
electric polarization P are induced by the mechanism proposed by Kimura et a/.- if the 
spin structure is canted relative to the c axis. This result offers a resolution to the is- 
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sue raised recently concerning the microscopic origin of the magnetolectric effect in this 
compound JL12i28«22 Such a canting occurs in crystals with rhomboderal symmetry only as a 
consequence of the trigonal interaction term F^. 

The method, and some of the results presented here, are also relevant to other frus- 
trated triangular antiferromagnets. Although hexagonal MnBr 2 shows planar anisotropy 
and strong inter-plane coupling responsible for a period-4 c-axis modulation, it shares the 
same basal-plane P4-IC transition as CuFeCv— NaFeC>2 has the same rhombohedral sym- 
metry as CuFeC>2 and a P4-IC transition but also with Sic and a more complicated P4 
spin structure.— Two of the three basal-plane triangular crystal directions show period-4 
modulations and the third is ferromagnetic. This can be described by q = (vr/2, n/A) xy or 
equivalently q = (rc/2, 7i/2) a b, which satisfy the umklapp requirement 4Q = G. Differences 
with CuFeC>2 could be due to additional inter-plane exchange interactions. It is of inter- 
est to note that this material does not appear to exhibit the structural phase transition to 
monoclinic symmetry found in CuFeC>2 below T/vi, providing another example where the 
stability of a P4 phase is not related to a structural distortion. Other potential applications 
of the present model include the hexagonal magnetoelectric antiferromagnet RbFe(Mo04)2 
which exhibits a complex phase diagram involving P3, P4 and P5 spin structures,— and a 
series of rhombohedral magnetoelectric antiferromagnets Ad02 (A=C r , A g , Li, or Na).— 
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